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Spacecraft  Sheath  AAodificatian  During 

Beam  Ejection 


1.  INTRODUCTION 

1  2 

In  recent  space  experiments,  *  the  potential  difference  created  between 
spacecraft  ground  and  the  ambient  plasma  during  the  ejection  of  a  beam  of  elec¬ 
trons  from  a  sounding  rocket  payload  in  the  ionosphere  has  been  found  to  be  much 

3 

less  than  had  originally  been  theoretically  predicted.  To  determine  the  reasons 
for  this  limited  potential  difference,  large  vacuum  chamber  tests  were  conducted 
in  which  electron  and  ion  currents  were  ejected  from  a  payload  into  a  simulated 
ionosphere. 
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Motivated  by  the  experimental  results,  sheath  ionization  models  '  for  small 
spacecraft  have  been  studied  for  a  plausible  explanation  of  the  observed  current- 
voltage  behavior.  When  an  electron  beam  is  emitted  from  a  spacecraft,  ambient 

g 

electrons  are  attracted  by  the  charged  spacecraft.  They  collide  with  the  neutral 
atmospheric  atoms  or  molecules  in  their  paths,  and  may  be  energetic  enough  to 

7 

ionize  the  neutrals  to  form  new  electrons  and  ions.  These  newly  created  charges 
alter  the  space  charge  current  arriving  at  the  spacecraft,  shifting  the  potential  to 
a  lower  value.  The  beam  electrons  are  assumed  to  be  energetic  enough  to  leave 
the  spacecraft  completely  and  to  play  a  negligible  role  in  the  ionization.  This 
mechanism  is  capable  of  explaining  the  low  potential  difference  observed,  and  a 
possible  non-monotonic  current-voltage  behavior. 


2.  MATHEMATICAL  FORMULATION 

The  method  of  approach  used  is  to  study  an  analytical  'plasma  probe' 

4  5  8  9 

model  •  '  *  with  space  charge  flow  of  electrons  accelerating  through  the  sheath 
surrounding  a  spherical  'probe'  that  represents  a  spacecraft  in  an  ionizable 
plasma  environment.  Magnetic  field  effect  is  ignored  in  this  model. 

The  beam  is  assumed  to  be  energetic  enough  to  leave  the  spacecraft  com¬ 
pletely,  and  is  not  stopped  by  its  own  space  charge  at  all.  As  the  beam  electrons 
leave,  the  spacecraft  becomes  charged  oppositely.  A  polarization  region  (sheath) 
is  formed  in  the  vicinity  of  the  spacecraft.  In  our  model,  ions  are  assumed  de¬ 
pleted  due  to  charge  repulsion  inside  the  sheath  (Figure  1). 

The  depletion  radius  r^  will  be  defined  by  the  balance  of  the  outgoing  beam 
current  with  the  incoming  ambient  current.  For  a  beam  current  1^^,  the  depletion 

radius  r  is  determined  by 
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Figure  1.  Sheath  Formation  During  Beam  Emission 
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where  is  the  thermal  velocity,  and  n^  is  ttie  number  density  of  ambient  elec¬ 
trons.  Some  typical  values  of  sheath  radius  as  calculated  by  means  of  Eq.  (1)  are 
shown  in  Figure  2. 
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Figure  2.  Parametric  Dependence  of  Sheath  Size 


The  potential  9  siny  point  inside  the  sheath  is  governed  by  the  Poisson 
equation 

*  -  -  f  . 

o 

where  p  is  the  space  charge  density,  and  is  the  permittivity  of  empty  space. 


.3.  SPHERICAL  SYMMETRIC  SYSTEM 


To  simplify  the  geometry,  we  assume  spherical  symmetry  in  the  spacecraft 
and  sheath  system.  Equation  (2)  becomes  simply  a  radial  equation: 


1  a_  /  2^ 

r^  3r  \  3r 


p(r) 

e„ 


where  the  gradient  of  the  potential  ^  gives  the  electric  field  E 


=  -  E(r)  . 


Taking  into  account  the  electron  and  ion  pairs  created  as  a  result  of  ionization,  the 
charge  density  p  at  any  point  r  in  the  sheath  is  given  by  the  sum  of  charge  densities 
(Figure  3): 

p(r)  =  e(n''’(r)  -  n'(r)  -  n^(r)]  ,  (5) 

where  n^  is  the  return  current  (primary)  electron  density,  while  n"*^  and  n  are  the 
ionization  ion  and  electron  densities,  respectively,  due  to  return  current  electron 
collisions  with  neutrals. 

The  ionization  electron  density  n  (r)  is  due  to  all  ionizations  that  occur  out¬ 
wards  of  r,  and  the  density  n^(r)  of  ions  at  r  is  due  to  all  ionizations  that  occur 
inwards  of  r.  Thus,  for  a  spacecraft  of  radius  R, 


[2e  U(r)  -  (i)(r')|/m„) 
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Figure  3.  Ionization  Pair-Creation  in  Sheath 


where  the  rate  of  ionization  ia  a  function  of  the  mean  free  path  X,  the  electron 
current,  and  the  velocity-dependent  probability  P  of  ionization: 


A’^P{Vg(r’)J  ng(r')Vg(r») 


(f 


Figure  4  shows  the  ionization  cross-section  of  atomic  oxygen  as  a  function  of 
energy.  It  has  a  peak  at  ~  90  eV  and  a  cutoff  at  ~  14  eV  for  typical  atmospheric 
neutrals. 


4.  NUMERICAL  METHOD 

To  solve  the  system  of  equations  [Eqs.  (3)  to  (8)],  one  divides  the  space  of 
the  sheath  into  N  concentric  shells,  and  sets  up  N  equations  for  the  N  unknowns 
0.  (see  Figure  5).  In  view  of  the  complexity  of  the  ionization  terms  in  Eqs.  (6) 
and  (7),  it  is  impossible  to  solve  these  equations  exactly.  Instead,  one  seeks  the 
approximate  solutions  that  minimize  a  function  F,  the  mean  square  of  f^^,  con¬ 
structed  from  the  radial  Poisson  equation  (Eq.  (3))  for  the  i-th  cell,  where  i=  1, 
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Figure  4.  Ionization  Cross-Section  as  a  Function  of  Electron  Energy 

f.(Ej,...EN)  =  (r^E).^^  -  (r^E).  -^(r2p(Ej....Ej^)].^r  . 

o 

where  the  electric  field  E  [Eq.  (4))  is  constructed  in  a  finite  difference  scheme 
<*.  =  Ar(E.+2E.^j+E.^2)/4  .  C 

The  numerical  method  used  to  solve  Eqs.  (8)  to  (10)  is  the  standard  Newton 
Ralphson  method  of  iteration: 


SPACECRAFT  POTENTIAL  (Voltt) 


To  start,  a  set  of  trial  solutions  is  used  in  the  Newton-Ralphson  iteration 
process,  and  a  convergent  set  of  solutions  is  sought  for  each  set  of  input  param¬ 
eters  such  as  beam  current,  ambient  electron  density,  ambient  electron  temper¬ 
ature,  mean  free  path,  and  spacecraft  radius. 


5.  RESULTS  AND  DISCUSSION 

Figure  6  shows  the  computed  results  of  spacecraft  potential  as  a  function  of 
electron  beam  current  for  various  electron  density,  electron  temperature,  and 
mean  free  path.  The  non-monotonlc  behavior  of  potential  current  curves  shows 
up.  At  low  currents,  the  potential  increases  with  beam  current.  When  the  cur¬ 
rent  increases  further,  ionization  occurs  inside  the  sheath.  The  potential  then 
turns  around  as  the  current  of  the  electron  beam  increases. 


Figure  6.  Typical  Non -Monotonic  Current  Voltage  Behaviors  for  (a)  Spherical 
and  (b)  Cylindrical  Spacecraft 

The  ion  and  electron  charges  created  by  ionization  alter  the  behavior  of  the 
space  charge  flow,  originally  governed  by  the  single  charged  Poisson  equation. 
The  potential  turns  to  a  lower  value,  and  stays  approximately  constant  as  current 
further  increases. 

In  the  rest  of  the  paper,  spherical  symmetry  will  be  assumed. 
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To  examine  the  mechanism  of  space  charge  modification  in  the  sheath  in 
greater  detail,  the  sheath  potential  profiles  for  various  beam  currents  are  dis¬ 
played  in  Figure  7.  As  the  potential  in  the  I-V  curve  turns  to  a  lower  value,  the 
potential  profile,  as  a  function  of  radial  distance,  shows  a  locally  flat  gradient. 
This  is  due  to  ions  created  inside  the  sheath  being  unable  to  move  out  quickly  due 
to  their  heavy  masses.  If  a  local  ion  charge  build-up  forms  a  potential  hump, 
ion  motion  would  be  two-way,  and  the  theory  would  then  break  down.  To  over¬ 
come  this  difficulty,  a  sweep  velocity  v^  is  added  to  the  ions,  and  Eq,  (7)  becomes 


n^(r) 


=  7/: 


fdnl 

UOr' 


•2  dr' 


[2e  U(r)  -  *(r')l/mj  +  v^  1 


27172 


(12) 
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Figure  7.  Relation  Between  Potential  Profile  and  I-V  Behavior 


It  is  argued  that  the  motion  of  a  spacecraft  relative  to  its  plasma  environment 
can  provide  such  a  sweep  velocity  v^  [Eq.  (12»  .  The  value  of  v^  is  of  the  order 
of  spacecraft  velocity,  and  is  an  arbitrary  input  to  the  computation.  However,  at 
a  higher  current,  a  potential  hump  again  shows  up,  the  computation  fails  to  con¬ 
verge,  and  the  technique  breaks  down.  It  is  conjectured  that  two-way  space 
charge  flows  should  be  accommodated  when  a  potential  hump  appears. 


The  effect  of  mean  ionization  distance  X  on  the  potential  of  a  spacecraft  emit¬ 
ting  an  electron  beam  is  shown  in  Figure  8.  In  the  region  where  ionization  is 

Q 

negligible  (X  «  lo  ),  the  potential  curves  are  flat,  and  the  potential  increases  witii 
beam  current.  As  the  mean  ionization  distance  decreases  (i.  e. ,  the  probability 
of  ionization  increases),  the  potential  curve  turns  down,  and  stays  approximately 
constant  as  the  ionization  distance  decreases  further. 


MEAN  IONIZATION  LENGTH  (m) 

Figure  8.  Spacecraft  Potential  as  a  Function  of  Mean  Ionization  Length 

Ambient  electron  density,  n^,  is  an  important  parameter  affecting  spacecraft 
potential.  In  Figure  9,  starting  with  large  values  of  electron  density,  the  poten¬ 
tial  increases  as  n^  decreases.  This  is  well  expected,  since  the  spacecraft  po¬ 
tential  depends  on  how  much  the  ambient  electrons  can  provide  in  the  incoming 
current.  However,  as  n„  decreases,  the  depletion  radius  r  increases  [r  oc  n'^^^ 
in  Eq.  (1)] ,  Hence,  the  total  distance  travelled  by  an  incoming  electron  increases, 
thereby  producing  more  ionization.  As  a  result,  the  potential  is  eventually  low  - 
ered  when  ionization  becomes  important. 

The  ambient  electron  temperature,  T  ,  plays  a  role  analogous  to  that  of  elec¬ 
tron  density,  n^,  but  in  an  opposite  manner.  At  high  temperatures,  the  ambient 
electron  current  is  high,  and  keeps  the  spacecraft  potential  low.  As  temperature 
decreases,  the  ambient  current  decreases,  but  the  depletion  radius,  r^,  increases. 
Then,  ionization  becomes  more  probable.  Eventually,  at  sufficiently  low  temper¬ 
ature  Tg,  the  potential  carve  turns  down  because  of  ionization  (Figure  10). 


PLASMA  DENSITY  (/m*) 

Figure  9.  Spacecraft  Potential  as  a  Function  of  Plasma  Density 

For  increasing  spacecraft  radii,  the  non-monotonic  current-voltage  behavior 
stiil  persists  (Figure  11).  However,  increased  spacecraft  radius  lowers  the  max 
imum  spacecraft  potential  induced  by  beam  emission.  Also,  the  amplitude  of  the 
difference  between  the  maiximum  potential  and  the  minimum  (beyond  the  turn¬ 
around)  diminishes.  Figure  12  shows  a  plot  of  the  envelope  of  maximum  and 
minimum  potentials  for  various  spacecraft  radii. 
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Figure  10.  Spacecraft  Potential  as  a  Function  of  Electron  Temperature 


For  a  given  beam  current  Ij^  [Eq.  (1)) .  the  sheath  surface  area  remains  con¬ 
stant  and  is  unaffected  by  the  increase  in  spacecraft  radius.  The  sheath  thickness 
(defined  as  the  sheath  radius  minus  the  spacecraft  radius),  however,  diminishes. 
As  a  result,  a  lower  spacecraft  potential  is  sufficient  to  attract  ambient  electrons 
through  the  sheath  for  the  compensation  of  electron  beam  current  leaving  the 
spacecraft. 

Beyond  the  turn-around  point  in  a  current-potential  curve,  the  minimum  po¬ 
tential  is  limited  by  the  minimum  energy  required  to  ionize  a  neutral  molecule  in 
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the  atmosphere.  Since  such  a  minimum  energy  is  generally  of  the  order  of 
14  eV,  the  minimum  potential  in  a  current -potential  curve  is  expected  to  approach 
14  eV  asymptotically,  depending  on  the  model  of  ionization  used.  For  the  same 
reason,  if  the  maximum  potential  induced  by  beam  emissions  is  below  ~  14  eV,  no 
non-monotonic  behavior  is  expected. 

Figure  12  shows  the  calculated  envelopes  of  the  maximum  and  minimum 
(beyond  turn-around)  potentials  for  various  spacecraft  radii  in  a  given  ambient 
environment.  The  amount  of  ionization  becomes  very  small  as  the  sheath  poten¬ 
tial  approaches  the  minimum  ionization  potential.  The  amplitude  of  the  potential 
drop  beyond  the  turn-around  also  approaches  the  value  of  minimum  ionization 
energy. 

There  is  another  critical  beam  current  that  manifests  itself  for  large  space¬ 
craft  but  not  for  small  ones.  This  current  is  determined  by  equating  the  sheath 
radius  to  the  spacecraft  radius.  If  the  sheath  radius  is  too  small,  the  spacecraft 
would  receive  enough  ambient  electrons  to  compensate  beam  emissions  without 
being  charged  up.  The  potential  of  the  spacecraft  in  this  case  is  that  of  natural 
charging.  Beyond  this  critical  current,  the  beam  emission  is  able  to  swing  the 
spacecraft  to  an  opposite  potential  and,  hence,  control  the  charging  of  the  space¬ 
craft.  This  phenomenon  shows  up  in  the  calculations  (Figure  11). 

In  the  model  studied,  as  the  radius  of  a  spacecraft  increases,  three  regimes 
of  physical  behavior  can  be  identified.  Figure  13  shows  these  regimes  clearly. 

The  potential  vs  spacecraft  radius  cm-ve  is  relatively  flat  in  the  small  radius 
regime.  This  is  the  regime  in  which  saturated  ionization  occurs,  i.  e.,  the  regime 
beyond  the  minimum  potential  in  a  current-voltage  curve.  The  second  regime  is 
characterized  by  the  presence  of  the  potential  maximum,  which  is  the  main  feature 
of  non-monotonic  behavior.  The  third  regime  occurs  when  the  spacecraft  is  so 
large  that  its  radius  exceeds  the  sheath  radius  (measured  from  th’’  spacecraft 
center)  for  a  given  current.  The  beam  loses  its  control  of  the  spacecraft  potential, 
and  natural  charging  dominates. 

For  ion  beam  emissions  from  a  spacecraft,  the  return  current  of  ions  is  not 
capable  of  turning  around  the  spacecraft  potential  This  is  because  the  cross- 
sections  of  ion  impact  ionization  of  neutral  atoms  or  molecules  are  much  coo  low, 
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Figure  13.  Non-Monotonic  Behavior  of  Spacecraft  Potential 
as  a  Function  of  Spacecraft  Radius,  for  a  Given  Electron 
Beam  Current 

in  contrast  to  the  electron  case.  The  current-voltage  characteristics  for  ion 

2 

beam  emissions  are  expected  to  be  monotonic,  as  demonstrated  in  space.  The: 
seems  to  be  no  bound  for  how  high  the  spacecraft  potential  can  be  induced  if  the 
current  of  the  ion  beam,  which  is  assumed  to  be  energetic  enough  to  leave  com 
pletely,  is  increased. 
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